Suppression of the structural phase transition and lattice softening 
in slightly underdoped Bai_^K^Fe2As2 with electronic phase separation 
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We present x-ray powder diffraction (XRPD) and neutron diffraction measurements on the slightly under- 
doped iron pnictide superconductor Bai_^K^Fe2As2, = 32 K. Below the magnetic transition temperature 
= 70 K, both techniques show an additional broadening of the nuclear Bragg peaks, suggesting a weak 
structural phase transition. However, macroscopically the system does not break its tetragonal symmetry 
down to 15 K. Instead, XRPD patterns at low temperature reveal an increase of the anisotropic microstrain 
proportionally in all directions. We associate this effect with the electronic phase separation, previously 
observed in the same material, and with the effect of lattice softening below the magnetic phase transition. 
We employ density functional theory to evaluate the distribution of atomic positions in the presence of 
dopant atoms both in the normal and magnetic states, and to quantify the lattice softening, showing that it 
can account for a major part of the observed increase of the microstrain. 

PACS numbers: 74.70.-b 25.40.Dn 61.05.cp 71.15.Mb 



I. INTRODUCTION 

The recent discovery of superconductivity (SC) in layered 
iron arsenides^'^-'^''*'^'^'^ served as a powerful impetus in the 
search for novel superconductors with high critical temper- 
atures. Within this new family of compounds, the record 
holders for the highest known are the electron-doped 
1111-compounds Gd^.^Th^FeAsO** and Srj.^Sm^FeAsF^ 
(both with optimal = 56 K), whereas among the so- 
called 122-compounds the highest of 38 K was reached 
in the hole-doped Ba^.^K^FcjAsj near the optimal doping 
of X ^ 0.5.^°'^^ The parent compounds (x = 0) of both 
types of arsenides order antiferromagnetically (AFM) below 
a spin-density-wave (SDW) transition at temperatures 
in the range between 140 and 200 K, as seen by neutron 
scattering^^'^'^'^'*'^^-^^'^'''^®'^' and local-probe methods, such 
as iiSR^°'^^'^^'^^ and ^'^Fe Mossbauer spectroscopy.^^-^'^ At 
low doping levels, this SDW transition is always accom- 
panied by a structural phase transition at > from a 
high-temperature tetragonal (T) to a low-temperature or- 
thorhombic or monoclinic structure, which manifests itself 
as a longitudinal splitting of the in-plane nuclear Bragg 
peaks (hhO) j both in neutron scattering and x-ray diffrac- 
tion experiments.12,13.14,15,16,17.18.23,24,25,26.27 Up to now, to 

the best of our knowledge, there have been no reports of 
any iron pnictides, neither among parent nor doped com- 
pounds, where magnetic order would be observed without 
the development of a structural distortion. 

Though it is commonly acknowledged that the magnetic 
and structural order parameters in iron pnictides are inti- 
mately coupled, the details of the relationship between the 
two phase transitions still remain a puzzle. On the one hand, 
in 1111-compounds (but not in 122-compounds)^^ the struc- 



gesting itself as the driving force for the magnetic anisotropy 
of the SDW phase. On the other hand, the experimentally 
observed structural distortion cannot be reproduced in non- 
magnetic calculations. Therefore most theories consider 
the SDW instability an intrinsic property of the electronic 
system, driven either by the nesting of the electron- and 
hole-like Fermi surface sheets'^"'''^'^^'^''''^'*'''^ or by the local 
superexchange interactions in the framework of the Heisen- 
berg model.^^'^^'^^'^'^''^° Both scenarios imply that the struc- 
tural phase transition occurs as a consequence of the AFM 
ordering, and its somewhat higher transition temperature 
is explained as a response to anisotropic AFM fluctuations 
that persist even above 7^,28,39,40,41 

In this paper, we combine neutron scattering and x-ray 
powder diffraction (XRPD) experiments, along with theoreti- 



1.5 



1.4 



1.3 



1.2 



I J * (72 V2 3)j magnetic Bragg intensity 

[ -^^ • • (1 1 0)t longitudinal width 



i i 



T \ i 



160 : 



120 ; 



80 



CD 

40 W 

CD 

o 







20 40 60 80 100 120 140 
Temperature (K) 

Fig. 1 (color online). Neutron diffraction data showing the 
temperature-dependent longitudinal broadening of the (1 1 0)^ 
nuclear Bragg reflection (•) overlayed with the intensity of the 



tural phase transition precedes the magnetic one, ' sug- (H magnetic Bragg peak (*). 



cal calculations, to study the interplay between the magnetic 
and structural phase transitions in a slightly underdoped 
122-compound Ba^.^K^FejASj (BKFA), = 32 K, in which 
the onset of a phase-separated magnetic order occurs at 
T^j = 70 K according to our recent study performed on the 
same samples.'*^ Our experimental evidence indicates that 
macroscopically the sample preserves its tetragonal symme- 
try down to 15 K, well below T^. Instead of the structural 
transition to an orthorhombic phase at low temperatures, 
seen in more underdoped BKFA samples, here the lat- 
tice reacts to the magnetic order only microscopically, by 
an increase of the microstrain as observed in our XRPD 
measurements, without a macroscopic breakdown of the 
lattice symmetry. We argue that such an effect is most prob- 
ably related to a softening of the lattice below the magnetic 
phase transition in comparison to the high-temperature non- 
magnetic state, whereas the phase-separated coexistence of 
twinned magnetic domains and the non-magnetic phase"*^ 
suppresses the structural phase transition beyond the exper- 
imentally detectable limit, in spite of a relatively high SDW 
transition temperature. 

II. SAMPLE PREPARATION 

The single crystals of BKFA used for the present study 
were grown using Sn as flux in a zirconia crucible sealed in 
a quartz ampoule filled with Ar. A mixture of Ba, K, Fe, As, 
and Sn in a wt. ratio of BKFA:Sn = 1:85 was heated in a 
box furnace up to 850°C and kept constant for 2-4 hours to 
soak the sample in a homogeneous melt. The cooling rate 
of 3°C/h was then apphed to decrease the temperature to 
550°C, and the grown crystals were then decanted from the 
flux."*^ Sample characterization by resistivity and dc suscep- 
tibility measurements'*^ revealed a sharp SC transition at 
^c, onset = (32 ± 1) K, reproducible among different samples 
from the same batch. The same samples have been exten- 
sively studied by angle-resolved photoelectron spectroscopy 
(ARPES)'*^''*'*''*'^ and muon-spin rotation (/xSR) ."^^'^^ 

The neutron diffraction measurements were done on a 
~ 30 mg single crystal with in-plane (hhO) and out-of-plane 
(OOZ) mosaicities better than 1.5° and 2.5°, respectively, as 
determined from the full width at half maximum (FWHM) 
of the rocking curves. A few smaller single crystals from 
the same batch were ground into powder for XRPD analysis. 
The sample was then prepared by sprinkling a small amount 
of the powder onto a flat brass sample holder. 

III. NEUTRON DIFFRACTION 

We have measured the longitudinal width of the (1 1 0)i- 
nuclear Bragg reflection as a function of temperature, which 
is plotted in Fig. 1 together with the intensity of the (| 1 3)^ 
magnetic Bragg peak. One can clearly see the broadening of 
the nuclear Bragg peak at low temperatures, with an onset 
at Tju, which perfectly follows the magnetic intensity, and 
amounts to ~ 20% as T ^ 0. The most straightforward ex- 
planation for such broadening would be a weak orthorhom- 
bic distortion that leads to a splitting of the peak that is 




Fig. 2 (color online). Panels (a) and (b) present XRPD data mea- 
sured at 300 K and 15 K, respectively, (i) Scattered x-ray intensity 
as a function of the diffraction angle 20 (A = 0.7 A) fitted to 
the tetragonal I4/mmm space group. For 20 > 17° the plots are 
enlarged by a factor of three. The fit includes a few wt. % of tetra- 
gonal ^-tin from the flux as an impurity phase and some reflections 
of the brass sample holder as indicated by the reflection markers 
in (ii). (iii) The difference A between the experimental points 
and the fitting curve. The insets show tensor surfaces representing 
the normalized anisotropic microstrain distribution along differ- 
ent crystallographic directions. The distance of the surface from 
the origin corresponds the squared full width at half maximum 
(FWHM) of the microstrain along the corresponding directions 
in real space. The x-z and x-y cross-sections of both surfaces are 
shown in panel (c) for comparison. 

masked by the experimental resolution, as was also previ- 
ously observed whenever the AFM order was suppressed 
either by doping, as in CeFeAsO0.94F0.06 low temperature 
[Ref. 13, Fig. 2(d)], or by temperature, as in the parent 
compound LaFeAsO at T = 138 K [Ref 12, Fig. 4 (inset)] . 
To check this interpretation, we have performed XRPD 
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measurements of the same samples, with subsequent ana- 
lysis of the microstrain anisotropy, which is known to be 
helpful in detecting minute structural distortions related to 
possible phase transitions.'*^''*^ 

IV X-RAY POWDER DIFFRACTION 

The XRPD data for the structure refinement were col- 
lected at room temperature and at 15 K, as shown in Fig. 2 
(a) and (b) . The sample was placed in a closed cycle cryo- 
stat. X-rays of 0.7 A wavelength were selected by a double 
Si(lll) monochromator. The wavelength and zero-point 
error were calibrated using 8 precisely measured peaks of 
the NBS1976 flat plate alumina standard. The diffracted 
beam was analyzed by reflection from a Ge(lll) crystal 
before a Nal scintillation detector Data were taken at each 
29 step of 0.005° from 3° to 38.6° at room temperature 
and 2° to 52° at 15 K. The sample was rocked during the 
measurement for better particle statistics. All data were 
normalized for storage ring current decay by an ionization 
chamber monitor. 

XRPD data were analyzed using the program TOPAS 
{Bruker-AXS} . Both high- and low-temperature data could 
be interpreted in terms of a tetragonal I4/mmm space group 
symmetry both at room temperature and at T = 15 K (see 
Fig. 2). As impurity phases, a few wt. % of tetragonal /3- 
tin from the flux and some reflections of the brass sample 
holder were included in the refinement. The analysis of the 
anisotropic peak broadening in the powder pattern due to a 
microstrain distribution was performed using the Cartesian 
parametrization by Leineweber."*^-"*® 

The lattice parameters of the sample, as deter- 
mined from XRPD by Rietveld refinement using the 
fundamental parameters approach of TOPAS, are 
a = b = 3.9111(1) A and c = 13.3392(6) A at room tempera- 
ture and a = b = 3.90075(7) A and c = 13.2476(3) A at 15 K, 
which corresponds to a 1.2% decrease in the unit cell vol- 
ume at low temperature. From the dependence of the lattice 
parameters on doping,^° the average potassium content of 
X = 0.4 could be determined, in agreement with the results 
of our energy-dispersive x-ray analysis. No evidence was 
found for an orthorhombic distortion of the tetragonal lat- 
tice at low temperature. This conclusion is based on the 
absence of any orthorhombic splitting of the Bragg reflec- 
tions and the refinement of the lattice parameters. The 
isotropic microstrain distribution in the (hfcO) plane also 
does not hint at an orthorhombic distortion. 

The microstrain distribution represents the statistics of 
the deviations Ad of the interplanar spacings from their 
average values, normalized by the average spacings d, i.e. 
of the strain e = Ad/d, over the investigated specimen as 
a function of the crystallographic direction. Tensor sur- 
faces representing the squared FWHM of the anisotropic 
microstrain distribution along different crystallographic 
directions are shown as insets in Fig. 2 (a) and (b), whereas 
panel (c) shows the x-z [tetragonal (ac) plane] and x-y 
[tetragonal (ab) plane] cross-sections of both surfaces. The 
largest microstrains of the crystalline lattice both at 300 K 



and at 15 K are found in the c-direction (\B^\^ = 0.9% and 
1.1%, respectively) as compared to the average in-plane val- 
ues of |BJ|| = 0.65% and 0.82%. The flowerlike shape of the 
x-z cross-section indicates a negative correlation between 
the in-plane (hfcO) and the out-of-plane (OOZ) directions, 
which agrees with the opposite changes of the a and c lat- 
tice constants upon the variation of doping.^" 

The low-temperature increase of the microstrain amounts 
to ~ 20% relative to the corresponding values at room tem- 
perature both in the c-direction and in-plane. In other 
words, to a good approximation the two tensor surfaces 
are geometrically similar to each other, which would not be 
expected in the case of a weak orthorhombic distortion, as it 
should instead broaden only the in-plane peaks. Moreover, 
at both temperatures no considerable in-plane anisotropy 
is observed [i.e. anisotropy in the x-y plane, see Fig. 2(c)], 
which would be a sign for the onset of an orthorhombic 
phase transition, e.g. for an incomplete orthorhombic reflec- 
tion splitting. Such an anisotropy (cf also Ref 48, Fig. 1) 
characteristically precedes tetragonal-to-orthorhombic struc- 
tural phase transitions in Pb304^* and La2Ni04.^^ 

This lets us conclude that the origin of the microstrain at 
both temperatures is not related to a macroscopic structural 
transition to orthorhombic symmetry, but rather should 
be attributed to an increase of the microscopic distortions 
of the lattice. The microstrain distribution quantitatively 
represents the response of the lattice to structural defects, 
such as chemical inhomogeneities or dislocations, which are 
unavoidable in any real material. Therefore an increase of 
the microstrain below the magnetic transition can either 
indicate that the lattice becomes softer, i.e. increases its 
response to the local stresses upon entering the AFM state. 




Fig. 3 (color online). The 2^/2ax2V2bxc supercell with 50% of 
the Ba atoms randomly substituted by the K dopants that we used 
in our density functional calculations. The calculated statistical 
distributions of the five interatomic distances, which are marked 
in the figure, are presented in Fig. 4. 
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Fig. 4 (color online). Histograms of the calculated interatomic dis- 
tances for the supercell as defined in Fig. 3. The plots at the top of 
each panel are a result of the normal-state (non-magnetic) calcula- 
tion, whereas the up-side-down plots below represent the magnetic 
ground state. The solid lines are fits to a normal distribution. 

or that the local stresses themselves increase, causing a 
proportional increase of the microstrain. 

In the studied compound, both mechanisms could be im- 
portant. On the one hand, in the case of lattice softening, 
one would expect its direct influence on the phonon mode 
frequencies. Indeed, such an effect has been reported in the 



phonon spectra of two similar 122-compounds: polycrys- 
talline Srg 5K0 4Fe2As2 and Cao.6NaQ4Fe2As2.^^ There, soft- 
ening of phonon modes below 10 meV has been observed 
by inelastic neutron scattering upon cooling from 300 K to 
140 K, despite the decrease of the unit cell volume at low 
temperature. More recently, softening and narrowing of 
several phonon modes below the spin density wave transi- 
tion was also observed by Raman scattering in underdoped 
Sri_^K^Fe2As2 and in the parent BaFe2As2 single crystals.^'* 
On the other hand, the phase-separated coexistence of AFM 
and paramagnetic phases in this material"*^ and the pres- 
ence of twin AFM domain boundaries^® should lead to an 
increase of local stresses below due to the magnetic 
anisotropy of individual AFM domains. 

To quantify the relative role of these two possible causes 
of the increased microstrain, we present here an estima- 
tion of the lattice softening across the magnetic transition 
based on our density functional calculations and show that 
it is comparable in magnitude and therefore could possi- 
bly provide a considerable contribution to the additional 
microstrain observed in the XRPD measurements. 

V DENSITY FUNCTIONAL CALCULATIONS 

Density functional calculations were performed using the 
projector augmented-wave^^'^^ method in the framework 
of the generalized gradient approximation. We have 
chosen a large 2V2a x 2V2b x c supercell, where 50% 
of the Ba atoms were substituted by K to model a random 
distribution of the dopants, as shown in Fig. 3. Using the 
Vienna ab initio simulation package (VASP) plane wave 
code,^''^°'^^-^^'^^ we have carried out the crystal structure 
optimization of the cell parameters and all ionic positions 
within the supercell to determine their displacements from 
the high-symmetry positions due to the introduced chem- 
ical disorder. The unit cell volume was fixed during the 
course of structure optimization to the experimental room- 
temperature value for the parent compound, as derived 
from a=b = 3.9625 A, c = 13.0168 A,^'* which also agrees 
well with our room-temperature value measured for BKFA. 
The cutoff energy of the plane-wave expansion was 367 eV, 
and the Brillouin zone sampling mesh was 4x4x4 with 
its origin at the F point. In the final optimized geometry, no 
forces on the atoms exceeded 0.01 eV/A. 

First, a non-magnetic calculation was performed, which 
represents the high-temperature (normal) state. As ex- 
pected, it revealed no deviations from the tetragonal symme- 
try, and resulted in the lattice parameters of a = b = 3.927 A 
and c = 13.258 A, which are reasonably consistent with the 
results of the XRPD structure refinement discussed above. 

The results of the crystal structure optimization are pre- 
sented in Fig. 4 (top plot in each panel), which shows his- 
tograms of the five interatomic distances, as defined in Fig. 3, 
fitted to a Gaussian distribution (solid lines). Naturally, the 
largest atomic displacements due to potassium substitution 
are observed in the Ba/K plane itself, where the distances 
^Ba/K between the neighboring Ba/K atoms vary by 0.176 A 
or 4.48% (0.0448 r. l.u., where r. l.u. stands for relative 
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lattice units), as estimated by the FWHM of the distribution. 
Within the FeAs block of layers, the out-of-plane atomic 
displacements (the buckling of the As and Fe planes) are 
the largest, and amount to 0.068 A (1.65% or 0.0051 r. l.u.) 
for the As layer and 0.088 A (1.32% or 0.0066 r. l.u.) for the 
Fe layer. In-plane distortions are notably smaller: 0.023 A 
(0.60%) and 0.013 A (0.45%), respectively 

Finally, we performed a spin-polarized calculation for 
the low-temperature striped AFM state. Collinear magnetic 
moments were self-consistently determined within the cal- 
culation. The corresponding histograms are shown in the 
same figure at the bottom of each panel. The most no- 
ticeable effect is the splitting of the Fe-Fe nearest neighbor 
distance, Xpg, which indicates the tendency of the system 
towards an orthorhombic distortion despite the presence of 
the dopants. As already mentioned above, such a transition 
is however suppressed macroscopically in the sample due 
to the presence of twin AFM domains and phase separation. 
In addition, one sees that the nearest-neighbor Fe-Fe inter- 
layer distance Zp^ increases, while that of As, z^s, decreases, 
which corresponds to the stretching of the Fe-As tetrahedra. 

Of more relevance for the present paper is the small but 
not negligible increase in the width of the distribution for 
every interatomic distance, as compared to the normal state, 
which we associate with the sought lattice softening effect. 
The variation of interatomic distances in the Ba/K plane 
increases by 3% to 0.182 A (4.67% or 0.0467 r. l.u.). The 
buckling of the As and Fe planes increases to 0.097 A (2.38% 
or 0.0073 r.l.u.) and 0.099 A (1.46%) or 0.0073 r.l.u.), 
respectively, which represents an increase by 43% and 12% 
relative to the corresponding normal-state values. In-plane 
distortions increase to 0.032 A (0.82%) for the As layer and 
0.019 A (0.67%) for the Fe layer (~ 40% of relative increase 
in the width). 

The observed changes in FWHM of the distributions 
between the AFM and the normal states are statistically 
significant and are observed consistently for all five con- 
sidered interatomic distances. For the FeAs block of lay- 
ers, they amount to ~0.002 r.l.u. on average both in- 
and out-of-plane. This is comparable with the increase 
of the microstrain |Bj^(15JC) - |Bj^(300ii:)= 0.002 and 
|B,|||(15JC) - |B,|||(300iif)= 0.0017 observed in our XRPD 
experiment. We can therefore conclude that the softening 
of the lattice associated with the SDW transition provides a 
major contribution to the observed effect. An additional con- 



tribution from the increase of the local stresses at the AFM 
domain boundaries due to weak local distortions within 
each domain can not be excluded, however. 



VI. SUMMARY AND DISCUSSION 

We have presented an example of an iron pnictide su- 
perconductor, which does not break its tetragonal crystal 
symmetry macroscopically upon entering the magnetically 
ordered state. This conclusion is based on x-ray powder 
diffraction measurements with subsequent analysis of the 
microstrain anisotropy. Instead, we have observed a low- 
temperature increase of the microstrain proportionally in all 
crystallographic directions, which has a magnetic origin and 
mostly originates from the softening of the crystal lattice 
below the SDW phase transition. A detailed analysis of the 
lattice structure in the presence of randomly distributed 
dopant atoms has been presented both in the normal and 
AFM states, confirming this conclusion. This does not ex- 
clude that a weak orthorhombic distortion possibly happens 
on a microscopic scale within each AFM domain, leading 
to an increase of the local stresses at the domain bound- 
aries, whereas the mesoscopic electronic phase separation"^^ 
suppresses the breakdown of the tetragonal symmetry on 
lateral scales larger than the typical size of the AFM domains. 
It can be argued that the observed magnetic state of the 
lattice represents a crossover between the well-developed 
orthorhombic phase at low doping levels and the normal 
tetragonal phase typical for the overdoped region of the 
phase diagram, where no magnetic transition occurs even 
at lowest temperatures. 
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